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ABSTRACT 

We investigate the weak gravitational lensing effect on the two-point correlation 
function of local maxima (hot spots) in the cosmic 21 cm fluctuation map. The intrinsic 
two-point function has a pronounced depression feature around the angular scale of 
6 ^ 40', which depends on the observed frequency and corresponds to the scale of 
the acoustic oscillation of cosmic plasma before the recombination. It is found that 
the weak lensing induces a large w-dependent smoothing at that scale where w is the 
equation of state parameter of dark energy and thus provides a useful constraints on 
the dark energy property combined with the depression angular scales on the two-point 
correlation function. 
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1 INTRODUCTION 

There has been a growing interest to use the cosmic 21 cm background fluctuations for a useful tool to study the so-called 
dark age of the universe. The emission and absorption lines of the 21 cm spin-flip transition of neutral hydrogen produce 
Cosmic Microwave Background (CMB) brightness temperature fluctuations of order ±10 mK. By scanning through redshift 
frequencies of 21 cm line, it is possible to observe the evolution of the neutral hydrogen density w ith time. Therefore, the 
observation of t he 21 cm fluctuations is expected as a promising probe of the reionisation history ijZaldarriaea et ah 2004; 
McQuinn et al In order to measure the fluctuations, there are planned low-frequency radio arrays, for example, 

Mileura Wide field Array (MWA), Low Frequency ArrajH (LOFAR) and Square Kilometer Arra}ll (SKA). Moreover, the 21 
cm cross-correlation with other observations, e.g, CMB, large scale structures and galaxies, also studied beca use the cross 



corre l ation with complementary observation gives more information t han their respective auto-correlations IjAlvarez et al 
20061 : iTashiro et al.l[20oi : lAdshead fc Furlanettolboosl : iLidz et allEoogl ). 



The 21 c m observations also have a potential to reveal the Universe at high redshifts before the reionisation epoch 
i Madau et al.l[l99 7: Tozzi ct al. 2000). The physics of the 21 cm flu ctuations at high redshif ts is understood better than at 
lower redshifts around the reionisation epoch (for a detailed review, see iLewis fc ChaUinojboOTi ). The 21 cm fluctuations at high 
redshifts trace the matter power spectrum and are calculated in the linear theory. The observations of the acoustic oscillations 
in the 21 cm fluc tuations are expected to be a probe of the composition and geometry of the Universe ( Barkana fc Loeb, 20051 : 
Mao fc Wulliooi ). 

In addition to these primordial fluctuations, there exist secondary 21 cm fluctu ations. Mo dification by weak gravitational 
lensing is considered as one of the main sources of secondary fiuctuations (Man del fc Zald arriaga 2006). The path of the 
emitted 21 cm photon is perturbed by weak gravitational lensing of large scale structures. Gravitational lensing modifies 
the primordial 21 cm fiuctuations as in the context of CMB. However, in contrast to the CMB case where there is only one 
source plane at the last scattering surface, the redshifted 21 cm fiuctuations provide excellent sources for gravitational lensing 
because of the existences of many different source planes. 



^ http:/ /www. haystack. mit.edu/array/MWA 

^ http://www.lofar.org 

^ http://www.skatelescope.org 
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In this paper, we investigate the weak lensing effect on the 21cm fluctuations. In orde r to reconstruct matter den- 
sity fluctuations along the 21 cm photon path^ some authors h ave been studied the effect ( Mandel fc Zaldarriaga 20061 : 
Zahn fc ZaldarriagalliooeljMetcalf fc Whit J2OO7I : ILu fc PenlboOSl ). For example, Mandel and Zeldarriga studied the effect on 



the angular power spectra as well as the intrinsic, three-dimensional power spectra of t he 21cm fluctuations during the era of 
reionisation and the effect on both spectra is less than a percent in the interesting scales ( Mandel fc Zaldarriaga 20061 ). Metcalf 
and White studied the lensed shear map power spectrum in 21cm fluctuat ions and pointed out th e potential of producing 
high resolution, high signal-to-noise images of the cosmic mass distribution dMetcaff fc WhitelboOTi ). 

Here we are interested in the two-point correlation function before reionisation. If adopting the Gaussian assumption 
for the primordial fluctuation fleld, it is known that the peak statistics can provide additional information about intrinsic 
distribution of hot spots that those pairs have some characteristic separation angles I Heavens fc Sheth 19991 ). In particular 
there is a pronounced depression feature in the two-point function around the angular scale of 6 ~ 40' depending on the 
observed frequency. The weak lensing then redistributes hot spots in the observed 21 cm fluctuation maps from the intrinsic 
distribution. We found that the effect has large influence in the depression feature, typically several percent contrary to the 
lensing effect on the power spectra. Moreover, the effect does not very much depend on the details of the nonlinear structure 
formation. Thus the detailed observation of the redshift dependence of the weak lensing effect as well as of the angular position 
in the two-point function offer a very useful mean of investigating dark energy equation of state. 

This paper is organised as follows. In Sec. II, we calculate and discuss the two-point correlation function of hot spots 
in the 21 cm line fluctuations. In Sec. HI, we give the formalisms the gravitational lensing effect on two-point correlation 
function. In Sec. IV, we present the result of lensed two-point correlation function and the effect of dark energy equation of 
state. Finally, we conclude in Sec. V. As the flducial cosmological model in this paper, we assume the ACDM cosmology with 
cosmological parameters, h — 0.7 (Ho — h x 100 km s~^ Mpc~^), f2,„ = 0.3, Q,a ~ 0.7 and erg — 0.8. 



2 TWO-POINT CORRELATION FUNCTION OF HOT SPOTS IN THE 21 CM LINE FLUCTUATIONS 



In the case of the CMB temperature anisotropy, two-point correlation functions of hot spo ts can be calculated fr om the 
angular power spectrum of the CMB temperature anisotropy under the Gaussian assumption ( Heavens fc Shet"hl 19991 ). Since 
it is a good assumption that the 21 cm fluctuations before the reionisation have a Gaussian statistic, we can calculate two-pint 
correlation functions of hot spots in the 21 cm fluctuations from their angular power spectrum in the same way as in the case 
of the CMB temperature anisotropy. 



2.1 Angular power spectrum of the 21 cm fluctuations 



The observed 21 cm fluctuations at A = 21(1 -f- Zobs) cm can be written as ([Zaldarriaga et al.ll2004[ ) 



r(n, Zohs) = T2l(Zobs 



driW2i{riohs - ?7)V'2i(n,T?), 



(1) 



where r\ is the conformal time, the subscripts obs and mean the values at the redshift z^hs and the present time, respec- 
tively, and W2\(r\o\is — r;) is a response function which characterises the bandwidth of an experiment and normalised as 
/^da;W2i(a;) = 1. In this paper, we assume that W2\(j]o\,s — f]) is the delta function for simplicity. In Eq. ([!}, T2\(^z) is a 
normalisation constant which is given by 



T2x{z) ~ 23 mK 



0.02 



) \\o.^h?) V 10 y 



1/2 



and ^21 is the dimensionless brightness temperature written as 

Tcmb (»?) 



'021 (n, 77) = a;H(n,r?)[l + 5b(n,T])\ 



1 



Ts(n,r,) 



(2) 



(3) 



where ^i, is the baryon density contrast, xh is the fraction of neutral hydrogen, Tg is the spin temperature of neutral hydrogen 
and Tcnib is the CMB temperature. For simplicity, we assume that the fraction of neutral hydrogen is a;H = 1 because the 
redshifts we are interested in are before the reionisation epoch. We also take another assumption that the gas temperature, 
Tg, is heated by stars or QSOs, and much higher than Tcmb- In this case, Tg is well coupled to Tg and Ts ^ Tcmb- As a result, 
Eq. ([Sj does not have the dependence on Tg. Under these assumptions, the 21 cm fluctuations are determined only by the 
baryon density contrast. 

Note that the second assumption, Ts ^ Tg, might be invalid at high redshifts, if there are not sufficient heat sources to 
make Tg enough high. In this case, Ts is set by the bal ance between Tg an d Tcmb, and the absolute value of the amplitude of 
the brightness temperature fluctuations depends on Ts (iMadau et al ] |l997l ). However, the absolute value of the amplitude does 
not affect the two-point correlation function in our analysis, since the threshold for a hot spot in this paper is deflned relative 
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to the dispersion of the fluctuations and the correlation function is determined by the gradient and the second derivative of 
the angular power spectrum as shown in iHeavens fc ShethI (|l999h . Therefore, our final conclusion is not significantly changed 
in the case without the second assumption. 

After expanding Eq. ((3]l in Fourier series and using Rayleigh's formula, we can obtain spherical harmonic coefficients of 
the 21 cm fluctuations, 



afm(^obs) = 4:Tv(-iY 



(2^ 



(1 + ffJ'^)SbkaY {k, Zohs)Y(^{k), 



where Yim{k) is a spherical harmonic function and aY{k,z) is a transfer function for the 21 cm line, 
af{k,z)=T2i{z)D{z)je[k{rjo~v)]- 



(4) 



(5) 



Here, je is the spherical Bessel function and D{z) is the linear growth factor of the density fluctuations. We included the factor 
(1 + //i^) in order to take into account of t he redshift-space distor tion by the bulk velocity flelds, which is called "Kaiser 
effect", with n = k n, and / = din D/d\na (jBharadwai fc Alill2004h . The angular power spectrum of the 21 cm fluctuations 
is obtained from 



C2l(2:obs,^) = {|a£m(Zobs)|^>- 



(6) 



2.2 Two-point correlation function 

For calculating the two-point correlation function of hot spots in the 21 cm line fluctuations, we deflne the number density 
fluctuations of hot spots in the 21 cm fluctuation map as 

Sn,,ie) = I^iiW^, (7) 
f^pk 

where ripk and are the number density and the mean number density of hot spots above a certain threshold p. The 
threshold v is written as = Ar2i/(T2i with the 21 cm fluctuations AT21 and their dispersion, 

a2i = (|Ar2i|^> = ^ J deC2iii). (8) 

In what follows, we set u — 1 as the threshold. 

The two-point correlation function of hot spots is the ensemble average of the number density fluctuations, 

Cpk-pk(e) = (5n(6>i)<5n(6»2)), (9) 

where \0i — 02\ = 0- The detailed calculation of the correlation function of hot spots from the angular spectrum is written in 



Heavens fc Sheth l|l999l) . Therefore, we do not give the detailed calculation and we only show some results here. 



The left panel in Fig. [T] shows the two-point correlation functions of hot spots in the 21 cm fluctuations. In this flgure, 
we set Zoha ~ 30. The baryonic oscillation brings the oscillatory feature between 20 and 50 arcmin in the correlation function. 
Since the 21 cm fluctuations do not have Silk damping, the correlation function does not damp below 10 arcmin unlike that 
of the CMB temperature anisotropy. 

Varying the observational wavelength, we can obtain two-point correlation functions at different 2obs- We plot the corre- 
lation function for different Zobs in the right panel of Fig. [T] When Zobs are varied from high to low redshifts, the position of 
the oscillation shifts to large angle. This is because the apparent angular diameter of the baryonic oscillation becomes large 
at low redshifts. In addition, decreasing the amplitude of the oscillation at low redshifts is also explained by the apparent 
angular diameter. The increment of the apparent angular diameter makes the baryonic oscillation smooth in the angular 
power spectrum of the 21 cm fluctuations. Since the two-point correlation function are related to the gradient and second 
derivative of the angular power spectrum, the smoothing of the baryonic oscillation in low 2;obs decreases the amplitude of the 
two-point correlation function. 



3 GRAVITATIONAL LENSING EFFECT ON TWO-POINT CORRELATION FUNCTIONS 

During traveling to an observer, 21 cm photons are deflected by the gravitational potential of the density fluctuations along 
the path as CMB photons. Therefore, we measure the lensed 21 cm fluctuations. This means that the number density of hot 
spots at the direction 9 on the observation map corresponds to that at 6* + 56* on the source plane, where 58 is the deflection 
angle generated by the weak gravitational lensing effect. Therefore, the observed number density fluctuations of hot spots, 
Sn'^^, are represented as 

Snil''{e) = Sn{z,e + se). (lo) 
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where 5n{z, 9) denote the number density fluctuations of hot spots at the redsh ift z. 
The dispersion of 56 for the angular separation 9 is written as ijSeliakll 19961 ) 



(11) 



where |0i — 02| = 6 and ( ) denotes the ensemble average. In Eq. (|lip . ctgl.o and crGL,2 are the isotropic and anisotropic 
contributions to the lensing dispersion, respectively, and given by 



o-GL,o(6') 



M 

2-K I T 



CglW[1- JoW], 



M. 



4L,2(e) = Y^GL{£)J2m, 

Here, Cgl(^) is the angular power spectrum of the deflection angle. 



xiv) 



(12) 



(13) 



where Pm is the matter density power spectrum and W(?7,r;obs) is represented as V[^(77, r^obs) = x('? ^ '?o6s)/x('?obs) where x is 
the comoving angular diameter distance. 



The effect of weak gravitational len sing on the correlation functions of hot spots was investigated by iTakada et al 



1 200(][ ) (see also Takada fc Futamasel 2001 ). The lensed two-point correlation function is expressed with the dispersions of the 
deflection angle as 



Cpk-pk(^) 



(5n°'==(0i)5n°'^=(02))|ei-e,|=e = (<5ri(0i + 50i)5n(02 + 502))|ei-e,|=e 
(27r)2 (27r)2^ lOn^^n^.^e ; 



2^' 



C, 



pk — pk 



1 - ^4L,o(e) ) Mm + ^^ol^^^{9)J2(l9) 



2 



where, in order to obtain the final expression, we used the Gaussian assumption of Jn^ , 

(5n,<5n,,) = (27r)'Cpk-pk(^)<5'(^ - ^'), 
and the following approximation, 

(e"-'''^-''^')|.i-e,|=e ^ 1 - ^ [<7SL,o(e) +cos(2<pO'^GL,2(e)] . 



(14) 



(15) 



(16) 



The angular spectrum of the unlensed correlation function of hot spots Cpk~pk in Eq. (|14p can be related with the unlensed 
correlation function Cpk-p"k ^''C^) as 



pk— pk^, 



(^) = 27r / d99Ql'^;r\Q)Mm- 



(17) 
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4 WEAK LENSING AND DARK ENERGY 



First, we calculate the weak gravitational lensing effect on two-point ©correlation functions in the fiducial cosmological model. 
The left panel in Fig. [2] shows the lensed and the unlensed two point correlation functions. The effect o f weak gravitational 
lensing is smoothing of t he oscillatory feature of two-point correlation function as in the CMB case (jTakada et al.ll2000l : 



Takada &: Futamasell200lh . Therefore, the effect of gravitational lensing arises prominently around 30-50 arcmin where there 
are the top and bottom of the oscillatory features. 

In this calculation, we have not taken into account the nonlinear evolution of matter dens ity fluctuations, although 



this gives the amplification of gravitational lensing on small scales in the CMB case l|Seliaklll99A . As shown in Fig. (2] the 



scale where the gravitational lensing effect arises prominently is about more than 30 arcmin. According to Eq. (|12p . most 
contribution of the gravitational lensing effect at each 6 comes from Cgl(^) where I corresponds to £ ^ n/O. Therefore, the 
gravitational lensing effect at about 30 arcmin is sensitive to Cgl(^) at ^ ~ 300. The discrepancy between the linear and 
nonlinear density power spe ctra on the power spectra of the deflection angle is made on small scales £ > 1000 (see Fig. 1 in 



Mandel fc Zaldarriaga 20061 ). Thus, neglecting the nonlinear evolution is valid in this paper. 



Next, we study the effect of the dark energy equation of state w on weak gravitational lensing in the two-point correlation 
function. Fig.|3]shows the lensed two-point correlation function for different w. One of the effects of w is the shift of the position 
of the oscillatory feature in the unlensed correlation function. Decreasing w means that the acceleration of the universe in 
the dark energy dominated epoch becomes high and the distance to ^obs increases. As a result, the oscillatory feature in the 
correlation function shifts to small angle scale as in the case of the baryonic oscillation in CMB or galaxy redshift surveys 
I Blake fc Glazebrookl [20o3 ) . Gravitational lensing does not affect the position of the oscillatory feature. Therefore, even in 



lensed correlation functions, the w-dependence of the position of the oscillatory feature is same as in the case of unlensed 
correlation functions. 

Varying the observational wavelength, we can obtain two-point correlation functions at different Zobs. The left panel of 
Fig. S] shows the evolution of the position of the oscillatory trough at the redshift ^obs for different w. As Zobs decreases, the 
apparent angle of the baryonic acoustic oscillation becomes large. Therefore, the position of the oscillatory trough shifts to 
large scale with Zobs decreasing. 

We present the evolution of the fractional change by gravitational lensing at the trough A^/^pJ^i^Tp^'^'^ with — (^pklpk — 
Cpk-pk"*)i a function of Zohs for different w in the right panel of Fig. U) The deflection angle depends on the distance to 
the redshift Zohs which is the redshift of a 'source plane', as described in Eq. (|13ll . Decreasing Zohs means that the distance 
becomes short and the deflection angle by gravitational lensing decreases. Therefore, the integrated value of Eq. (|13p is smaller 
for low Zohs than for high Zohs- 

The amplitude of A^/^^;,'!.™^'^ depends on w as shown in the right panel of Fig. |4l Decreasing w with keeping Zobs makes 
the comoving distance to Zobs large. As a result, integral range of Eq. (I13|l become large and the gravitational lensing effect is 
enhanced for small w. The modification of the growth rate of density fluctuations by the dark energy equation of state affects 
weak gravitational lensing. However we found that this effect is minor in our parameter region ( — 1.2 < w < —0.8), compared 
with the effect of the modification of the distance to Zohs- 

At the last, we investigate which redshift makes the most contribution to the lensing effect. Since the gravitational lensing 
effect mainly arises at about 30 arcmin in the 21 cm fluctuations between Zobs ~ 20 and 40, the dominant contribution to the 
lensing dispersion in Eq. (|12[l comes from Cgl{£) at f ~ 300. We plot the redshift distribution of Cgl(^) a,t £ = 300 in Fig.O 
Varying Zobs does not make the distribution change radically, because the radial distances to the source plane at each Zohs are 
not different strongly. For all Zobs, the distribution has a peak around z ~ 1.5. Therefore, gravitational lensing effect on the 
21 cm fluctuations is sensitive to the matter density fluctuations at this redshift 



5 CONCLUSION 

We have calculated two-point correlation functions of hot spots in 21 cm fluctuations and studied the weak gravitational 
lensing effect on the correlation function. Particularly, we have examined its possibility as a probe of the equation of state of 
dark energy w. 

The two-point correlation function of hot spots in the 21 cm fluctuations is more smoothed than that of the CMB 
temperature anisotropy. On large scales, the correlation function is very flat. However, the amplitude of the correlation 
function increases toward small scales due to the absence of Silk damping in the 21 cm fluctuations. The baryonic acoustic 
oscillation produces the oscillatory feature around 30-50 arcmin. By decreasing the observation redshift Zohs, the oscillatory 
feature shifts to large scales. Therefore, the position of the oscillatory feature in multi-frequency observations plays a important 
role in the decision of the cosmological parameter as the baryonic oscillation in the CMB temperature anisotropies and large 
scale structures. 

The weak gravitational lensing effect on the two-point correlation function appears on the oscillatory feature. The effect 
is smoothing of the feature by smearing the baryonic oscillation without shifting its position. The advantage of the 21 cm 
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30 40 50 60 70 SO 29 30 31 32 33 34 35 

9 (arcmin) 9 (arcmin) 

Figure 2. The two-point correlation functions of hot spots in the 21 cm fluctuations. In the left panel, the solid line represents the 
lensed correlation function and the dotted line indicates the unlensed correlation function. We set Zo\,s = 30 and erg = 0.8. The right 
panel shows the result around 30 arcmin in the same case. 




Figure 3. The lensed two-point correlation functions of hot spots for different w. The dotted, solid, and dashed lines represent the lensed 
correlation functions for w = —1.2, w = —1.0, and w = —0.8, respectively. For comparison, we also plot unlensed correlation function as 
the thin lines. We set ^o^s = 30 a^nd ""S = 0.8. 



fluctuation observation is that we can obtain independent lensed maps at different redshifts. As Zohs decreases, the distance to 
the source plane becomes small and the deflection angle by gravitational lensing decreases. As a result, the difference between 
the lensed and the unlensed correlation functions becomes small at low Zobs- 

We have studied the sensitivity of the 21 cm two point correlation function to the dark energy equation of state w. 
The effects of w on the two point correlation function appear on the shift of the position of the oscillatory features and the 
smoothing by gravitational lensing. Since the distance to Zobs depends on w, decreasing w makes the position shift to small 
scale and the smoothing emphasised. 




15 20 25 30 35 40 15 20 25 30 35 40 

"^obs '^obs 

Figure 4. {Left) The evolution of the position of the oscillatory trough for different w. The dotted, solid, and dashed lines are for 
w = —1.2, w = —1.0, and w = —0.8, respectively. {Right) The evolution of fractional changes of the two-point correlation function for 
different w. The dotted, solid, and dashed lines represent the evolutions for w = -1-2, w = —1.0, and w = —0.8, respectively. 
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Figure 5. The redshift distribution of Cgl{^) at ^ = 300 as a function of z^^^^. The dotted, sohd and dashed lines are for z^^^ = 40, 
■^obs = 30 and z^hs = 20, respectively. 



The evolution of the fractional change by the gravitational lensing effect, which is obtained from observations of different 
redshift slices, is useful for the constraint on w. The 21 cm fluctuations with gravitational effect can be estimated from the 
linear theory. Thus, it will be easy to compare these results with observational data and to obtain the constraint on w. 

We mention the comparison between the angular power spectrum and th e two-point correlation f unctio n of the 21 cm 
fluctuations for the detection of the gravitational lensing effect. According to Mandel fc Zaldarriaga ( 20061 ) . the fractional 
change by gravitational lensing in the angular power spectrum is about 1 %. On the contrary, the fractional change is 
enhanced to about 4 % in the correlation function. Therefore, the two-point correlation function is a better probe of the 
gravitational lensing effect than the angular power spectrum. 

Finally, we give some comments on observational aspects. In this paper, we focus on weak gravitational lensing on the 21 
cm fluctuations before the reionisation epoch (zobs > 15). Therefore, we neglect any effect of reionisation process on the 21 
cm fluctuations. Planned observational projects in the near future are aimed at the measurement of the 21 cm fluctuations 
during reionisation (zobs ~ 10 for LOFAR and Zoba ~ 15 for SKA). The 21 cm fluctuations from the reionisation epoch are 
studied with numerical si mulations by man y authors. The angular power spectrum of the 21 cm fluctuations depends on the 
reionisation process, (e.g. Baek et al. 20081 ). Therefore, in order to measure w by the weak gravitational lensing effect in the 
near future observations, we must take into account the reionisation process precisely with numerical simulations. 
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